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RESUMEN 

Presentamos los primeros resultados de un proyecto dedicado al estudio combinado de estrellas masivas y 
sus regiones H n asociadas en la Galaxia. Mediante el estudio de regiones ionizadas por una unica estrella 
prctcndemos comprobar la validez de las prediccioncs de la nueva generation de modelos de atmosfera de 
estrellas masivas en terminos de distribution de flujo ionizante. En nuestros analisis considcramos el efccto 
de las posiblcs distribuciones de densidad del material ionizado. Para ello hacemos uso de varios tipos de 
observaciones, tanto estelares como de la region Hn. 

ABSTRACT 

We present the first results of a project aimed at the combined study of massive stars and their surrounding 
nebulae by means of a detailed study of Galactic Hn regions ionized by only one massive star. With this, we 
intend to check the validity of the new generation of massive star model atmosphere codes in terms of ionizing 
flux distribution. We take into account the effect of the nebular density distribution in our analyses. Various 
types of stellar and nebular observations have been collected for this purpose. 
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1. INTRODUCTION 



STARS: INDIVIDUAL: HD 37061 — H II REGIONS 
ISM: INDIVIDUAL OBJECTS: M 43 



ISM: 



The intense far ultraviolet radiation emitted by 
early OB-type stars ionizes the interstellar medium, 
generating the so-called HII regions. These ionized 
regions can be used to derive properties of the as- 
sociated stellar population (e.g. star forming rates, 
age). However, since the properties of HII regions 
crucially depend on the spectral energy distribution 
(sed) of the massive star population, and this part 
of the stellar flux is generally unaccesible to direct 
observations, the predictions resulting from massive 
star atmosphere codes are a crucial ingredient. 

The new generation of nlte, line blanketed 
model atmosphere codes (Hubcny & Lanz, 1995, 
Hillier & Miller 1998, Pauldrach et al. 2001, Puis 
ct al. 2005), which include a more realistic descrip- 
tion of the physical processes characterizing the 
stellar atmosphere, produce quite different ionizing 
SEDs than the previous plane-parallel, nlte/lte, 
hydrostatic models (Mihalas & Auer 1970, Kurucz 
1992, Kunze 1994). Some notes on this, and on 
the consequences on the ionization structure of H II 
regions, can be found in Gabler et al. (1989), Rubin 
et al. (1995), Najarro et al. (1996), Scllmaier et al. 
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(1996) and Stasinska & Schaerer (1997). Although 
new predictions seem to go in the right direction 
(viz. Giveon et al. 2002, Morisset et al. 2004), 
non-negligible differences can be found between the 
different stellar atmosphere codes (viz. Martins et 
al. 2005, Puis ct al. 2005). 

We arc performing a study of Galactic Hn re- 
gions ionized by only one massive star to check 
the validity of the ionizing SEDs predicted by the 
new generation of massive star model atmosphere 
codes. In forthcoming papers we will present the 
complete combined study of the various nebulae and 
their ionizing stars; here we show some preliminary 
results from our study of M43. 



2. OBSERVATIONAL DATASET 

Stellar and nebular observations were carried out 
with the Wide Field Camera (wfc; narrow-band 
imaging in Ha, H/3, [O III] and [S II] ) and the Inter- 
mediate Resolution Spectrograph (ids; spectroscopy 
of the ionizing star in the range 4000 - 5000 A + the 
Ha region, and long slit nebular spectroscopy in the 
range 3600 -9700 A); both instruments attached to 
the Isaac Newton Telescope (int) . The long slit neb- 
ular observations were divided into smaller apertures 
along the nebular radii. 
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Fig. 1. Fitting of fastwind synthetic profiles (broadened to «sini = 2f0 
kms -1 ) to the optical H and He lines of HD 37061. Two sets of stellar 
parameters arc shown to illustrate the accuracy of the stellar parameters 
determination. 
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Fig. 2. Spatial variation of various nebular line ra- 
tios resulting from three CLOUDY constant density 
spherical models. Each model uses a different pre- 
dicted SED. FASTWIND models are the ones used in 
Figure 1 (same color code); a TLUSTY model with 
the same stellar parameters as the FASTWlNDmodcl 
which best fits the optical HHe lines (Figure 1) has 
also been considered for comparison. Nebular ob- 
servations from long slit spectra are also included. 



3. STELLAR AND NEBULAR ANALYSES 

M43 is an apparently spherical Hn region, ioni- 
zed by HD 37061 (classified as B1V, though our 
spectrum clearly shows the He II 4541 line, indicat- 
ing that the star is rather BO V-B0.5 V). The stellar 
parameters of HD 37061 were obtained by visual fit- 
ting of FASTWIND (Puis et al. 2005) synthetic profiles 
to the optical H and He lines (see Figure 1). Once 
the SED resulting from the fastwind model was re- 
scalcd to fit M v = 3.5 (from the stellar photometry 
and considering a distance of 450 pc), a value of 
log Q(H°) = 47.2 could be derived. This later value is 
in good agreement with the nebular H a luminosity 
calculated from the wfc H a image, indicating that 
the nebula is ionization bounded. From the nebu- 
lar [S II] 6731/6716 line ratio we inferred N c = 550 
cm . In a first approach, we have used the nebular 
abundances derived by Rodriguez (1999) 

4. PHOTOIONIZATION MODELS 

Figure 2 illustrates three of the diagrams we use 
to compare the predictions of photoionization mod- 
els (cloudy, Ferland et al. 1998) with the observed 
nebular constraints. The first one is an indicator of 
the nebular temperature, while the other two illus- 
trate the nebular ionization structure. These plots 
compare three constant density spherical CLOUDY 
models: two of them take the SEDs from FAST- 
WIND models with different stellar parameters, while 
the third one was calculated using the SED from a 
TLUSTY (Hubeny & Lanz, 1995) model with the same 
parameters as the FASTWIND model which best fits 



the optical HHe lines. The differences between the 
photoionization models arc important. Note, how- 
ever, that geometries other than spherical and non- 
constant density distributions can affect the result 
(Morisset et al. 2005). This will be also taken into 
account when inferring which model better fit the 
nebular constraints. In forthcoming papers we will 
present the complete analysis in detail. 
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